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Abstract 

A model of a lossy nonlinear fiber grating with a "hot spot", which combines a local gain and 
an attractive perturbation of the refractive index, is introduced. A family of exact solutions for 
pinned solitons is found in the absence of loss and gain. In the presence of the loss and localized 
gain, an instability threshold of the zero solution is found. If the loss and gain are small, it is 
predicted what soliton is selected by the energy-balance condition. Direct simulations demonstrate 
that only one pinned soliton is stable in the conservative model, and it is a semi — attractor: 
solitons with a larger energy relax to it via emission of radiation, while those with a smaller energy 
decay. The same is found for solitons trapped by a pair of repulsive inhomogeneities. In the model 
with the loss and gain, stable pinned pulses demonstrate persistent internal vibrations and emission 
of radiation. If these solitons are nearly stationary, the prediction based on the energy balance 
underestimates the necessary gain by 10 — 15% (due to radiation loss). If the loss and gain are 
larger, the intrinsic vibrations of the pinned soliton become chaotic. The local gain alone, without 
the attractive perturbation of the local refractive index, cannot maintain a stable pinned soliton. 
For collisions of moving solitons with the "hot spot" , passage and capture regimes are identified, 
the capture actually implying splitting of the soliton. 

PACS numbers: 
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I. INTRODUCTION 



Solitons in any physical medium are subject to attenuation due to dissipation, hence 
it is necessary to apply gain which can support the solitons. A possibility which may find 
interesting physical applications is to create a localized gain, which will be a trap for solitons 
in a lossy medium. In terms of optical solitons, this can be easily realized for spatial solitons 
in a planar waveguide, where the gain may be applied to a narrow strip. However, for spatial 
solitons gain is a redundancy, as they are supported simply by the energy fiux through the 
soliton itself. Besides that, in this work it will be shown that, in the simplest case when the 
spatial solitons are governed by the nonlinear Schrodinger (NLS) equation, a pinned soliton 
supported by the localized gain can never be stable. On the other hand, realization of the 
local-gain trap for usual temporal solitons is impossible, as such a soliton runs (for instance, 
in an optical fiber [||]) with the group velocity of light. 

An unique possibility to create a gain-induced trap in a lossy medium is offered by a fiber 
grating, i.e., a Bragg grating (BG) written on an optical fiber. Fiber gratings are a basis for 
many photonic devices 0]. A challenging possibility is to use fiber gratings for the creation 
of pulses of slow light, which is a topic of great current interest [Q]. The possibility of the 
existence of slow pulses suggest to try a local-gain-induced trap for solitons. 

In fiber gratings, solitons exist due to the interplay between the Bragg refiection and 
Kerr nonlinearity of the BG-carrying fiber [^]. These solitons were predicted analytically 
and then they were created in the experiment [0]. Except for the case in which 
the BG soliton is very broad this species of optical solitons is distinct from the usual 
nonlinear-Schrodinger (NLS) solitons in ordinary nonlinear optical fibers. 

As it was mentioned above, search for very slow solitons in fiber gratings, which were 
predicted long ago 0, is an issue of great interest 0. The standard mathematical model 
of the nonlinear fiber grating predicts a whole family of zero- velocity solitons 0, ^ a 
part of which is stable 0. If realized experimentally, such a soliton would represent a 
pulse of standing light, with its left- and right-traveling components being in a permanent 
dynamical equilibrium. 

The BG solitons which have been observed in the experiment up to date are moving ones, 
their velocity being ~ 75% of the maximum group velocity of light in the fiber [0 . The quest 
for zero- velocity solitons may be facilitated by means of a local defect in the BG that exerts 
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an attractive force on a soliton, having thus a potential to be a sohton trap (note that a 
local defect may trap light in the fiber grating via the four-wave mixing without formation 
of a soliton [|12|). Besides its profound physical importance |]13[, such a soliton trap is also 
promising for the fiber-sensing technology ||14|| . 

The interaction of the soliton with an attractive defect in the form of a local suppression 
of BG was studied recently in Refs. |]TE[ and |M. The latter work also considered a trap 
combining the local BG suppression and a change in the refractive index (these two effects 
may come together as a manifestation of a local inhomogeneity in the BG-carrying fiber). 
As a result, it has been demonstrated that a gap soliton may indeed be pinned by the local 
defect, and, moreover, a narrow (delta-functional) defect uniquely selects parameters of the 
stable trapped soliton [|14| (this feature is described in more detail below). 

However, a trapped soliton will be destroyed by fiber loss. Indeed, taking into regard 
that the best fiber gratings (used as dispersion compensators in optical telecommunications 



T5| ) have the attenuation rate ~ 0.2 dB/cm (in hybrid grating waveguides, using glass and 



a sol-gel material, the attenuation may be lowered to 0.1 dB/cm |T^), it is easy to estimate 
that a standing soliton will be destroyed during the time ^ 5 ns. Therefore, it is necessary to 
support the trapped soliton by means of a locally applied optical gain, which is tantamount 
to the above-mentioned soliton trap induced by local gain. A related problem is a possibility 
to capture a moving soliton by the local gain. These issues are the main subjects of the 
present work. Besides the fundamental interest to having permanently maintained trapped 
optical solitons, they may also be interesting for applications, such as optical memory. 

The local gain can be provided, for instance, by a short resonantly doped segment in 
the BG-carrying fiber. In this connection, we note that the moving BG soliton that was 
observed for the first time in the fiber grating had the temporal width ~ 200 ps [0. Taking 
into regard the Lorentzian contraction (as it was mentioned above, the soliton was moving 
at a velocity equal ~ 75% of the maximum group velocity), the same soliton, if stopped, 
would have the spatial width on the order of a few centimeters. Thus, the locally pumped 
segment of the fiber may be approximated by a delta-function (which we assume in the 
analytical model below) if, roughly speaking, its size is ~ 1 mm. Note that the maximum 
gain that can be provided by an Er-doped amplifier is ~ 5 dB/cm |T^. Comparing this 
to the above-mentioned minimum damping rate 0.1 dB/cm, we conclude that the proposed 
scheme may be self-consistent. 
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Thus, the model to be formulated below assumes the local gain in the form of the delta- 
function (in the numerical part of the work, a smooth approximation to the delta-function 
is used). The model also includes an imaginary part of the localized gain, which actually 
accounts for a local perturbation of the fiber's refractive index (note that the dopant, if used 
to induce the local gain, may indeed affect the local refractive index). The imaginary part 
is necessary, as it will be seen that pure gain cannot maintain a pinned soliton in a stable 
state. As for the size of the refractive-index perturbation 5n, it will be seen that ^ra ^ 0.1 
is definitely sufficient to stably trap a soliton. In this respect, it is relevant to mention that 
a dopant added to a silica fiber usually induces a refractive index change 5n ~ 0.05, while 
in a polymer fiber it is ~ 0.2. The latter value can readily produce a stable pinned state of 
the soliton, and the former one may be sufficient too. 

A relatively large jump of the refractive index between the doped segment and the rest of 
the fiber may induce an additional effect, viz., reflection of light, although the reflection may 
be smothered by a sufficiently smooth profile of distribution of the dopant. An additional 
term in the basic model [see Eqs. (|l]) and (0) below] , induced by the refiection, would be the 
same which formally accounts for a local perturbation of the BG refiectivity. The latter type 
of the conservative perturbation was considered in Ref. fl^)) a conclusion being that its 
effect is quite similar to that directly produced by the refractive-index perturbation, which is 
directly included in Eqs. (|I|) and (^. For this reason, we do not consider the local reflection 
as a separate perturbation in this work (in any case, it can be easily added, if necessary). 

The paper is organized as follows. In section 2 we formulate the model, and produce some 
analytical results. First of all, we consider the stability of the zero solution against small 
perturbations in the presence of the uniform loss and localized loss. An instability-onset 
threshold is found, and it is demonstrated that the instability does not take place unless the 
gain has an imaginary part. For the comparison's sake, we also briefly consider an allied 
model problem, viz., the instability induced by a "hot spot" in the lossy NLS equation, which 
produces quite similar results. Then, we produce a family of exact analytical solutions for 
solitons pinned by the local inhomogeneity of the refractive index, in the absence of loss and 
gain. Another analytical result predicts what particular pinned soliton is selected by the 
balance of the loss and gain, provided that they are small perturbations. A soliton selected 
by the energy balance is also found in the lossy NLS equation with the hot spot (in that 
case, it is always unstable). 
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Section 3 reports results of direct simulations of pinned solitons, in both the conservative 
version of the model, and in the full one. A surprising result is that, in the conservative 
model, a single pinned soliton, out of their continuous family, is stable. If the initially 
created soliton has larger energy, it relaxes, via emission of radiation, to the stable one, and 
if the energy of the initial soliton is smaller than that of the stable one, the pulse decays 
into radiation; the former observation can be explained by known results for the stability 
of the ordinary solitons in the unperturbed BG model ||10|, A similar result is also 



obtained, in the framework of the conservative model, at the end of section 3 for a pair of 
repulsive defects: in the case when they can hold a soliton at the midpoint between them, 
any pulse either relaxes to a uniquely selected soliton, or decays (if the pair of repulsive 
inhomogeneities cannot hold a soliton, it splits the soliton in two). It is relevant to mention 
that a nonlinear BG structure with two local defects was very recently studied as a promising 
model for all-optical switching . 



In the full model, all the pinned solitons demonstrate persistent intrinsic vibrations; 
depending on the ratio of the loss and gain, and on the strength of the attractive refractive- 
index inhomogeneity, the amplitude of the vibrations may be small or large. In the former 
case, the pinned soliton may be regarded as an approximately stationary one, then the above- 
mentioned analytical prediction of the soliton selected by the balance between the loss and 
gain gives an error between 10% and 15%, which may be explained by extra radiation loss. 

In section 4, collision of a moving soliton with the hot spot is considered by means of 
direct simulations. As a result, regions of passage and capture are identified in the soliton's 
parameter space in both the conservative and full models. In fact, the capture is incomplete: 
a part of the soliton's energy gets trapped, giving rise to a pinned soliton, while the other 
part passes and rearranges itself into another soliton. If the inhomogeneity is strong, a 
conspicuous part of the energy may bounce back. 

II. THE MODEL AND ANALYTICAL RESULTS 
A. The model equations 

A localized gain inserted into a fiber grating is modelled by a 5 - function term added to 
the standard BG model, which includes coupled equations for the amplitudes of the right- 
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and left-travelling electromagnetic waves, u{x,t) and v{x,t) Q]. The full model is then 



u 



-i'yu + iT6{x) ■ u, 
-i'-yv + iT6{x) ■ v, 



(1) 
(2) 



iut + iux + f + {\u\^ + 2|t>| 
ivt — ivx + u + (^{vl"^ + 2\u\ 

where the maximum group velocity of light is normalized to be 1, the nonlinear terms 
account for the self- and cross-phase modulation induced by the Kerr effect, the linear 
couplings represent the mutual conversion of the waves due to the resonant Bragg scattering 
(the conversion coefficient is also normalized to be 1). On the right-hand side of Eqs. (|lD 
and (0), the fiber loss parameter 7 is real and positive, while the local-gain strength F may 
be complex, 

r = Ti + zFa , (3) 

its positive real part being the gain proper, while the imaginary part accounts for a localized 
perturbation of the refractive index (note that r2 > corresponds to a local increase of the 
refractive index). 

A stationary solution, which represents a soliton trapped by the "hot spot" , is sought for 

as 

u{x, t) = U (x) exp {—it cos 6) , v{x, t) = V{x) exp {—it cos 9) , (4) 

where ^ is a parameter of the soliton family. The substitution of Eqs. (^ into Eqs. (|l|) and 
(0) leads to equations 



. d 

dx 
d 



+ cos^ + Z7 + iV5{x) + (|f/p + 2\V\- 



-i^ + cos^ + Z7 + iT6{x) + + 2\U\'^ 



U + V 

v + u 



0, 
0. 



(5) 
(6) 



B. Stability of the zero solution 



We start the analysis with the linearized version of Eqs. ([T|) and (^, in order to analyze 
the stability of the zero solution in the presence of the localized gain. An eigenmode of small 
perturbations is sough for as 



u{x, t) = exp{—ixt — Kx) , v{x, t) = exp{—ixt — kx) at x > 0, 
u{x,t) = A_exp{—ixt + Kx) , v{x,t) = B_ exp{—ixt + kx) at x < 0, 
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(7) 
(8) 



where the frequency x "^^Y be complex, its imaginary part being the instabihty growth rate, 
and K must have a positive real part. Substituting this into the linearized equations yields 

K = ^l-(;^ + Z7)2, (9) 

where the square root is defined so that its real part is positive, and B± = 
— [{x + ^7) =F ii^] Further, the integration of the linearized equations in an infinitesi- 
mal vicinity of x = yields = A_ e^, = B_ e~^. 

Eliminating the amplitudes A± and B± by means of these relations, we obtain an equation 

X+ll^ = e-(r^+^r^), (10) 
X + il + 

where the expression @ was substituted for V. Finally, combining Eqs. (|10D and @, we 
can find the eigenfrequency sought for: 

X = -il + sgn (sin ■ cosh (Fi + 1^2) , (11) 

where the sign multiplier in front of the second term provides for the fulfillment of the 
above-mentioned condition Re k > 0. 

A straightforward consequence of Eq. ([TT| ) is an expression for the instability growth rate. 



Imx = —7 -|- (sinh Fi) | sin F2I . (12) 

Thus, the localized gain gives rise to the instability of the trivial solution, provided that it 
is strong enough: 

sinhFi>sinh((FOJ = ^-^. (13) 

I sni i 2 1 

Note that the instability is impossible in the absence of the local refractive-index perturbation 
F2. The instability-onset condition (p!3| ) simplifies in the limiting case when both the loss 
and gain parameters are small (while F2 is not necessarily small): 

Fi>(Fi)^^^— (14) 
sm 1 2 



Lastly, we notice that the instability is oscillatory: as it follows from Eq. (|TT|), Rex 7^ 0, 
unless cos F2 = 0. 

As the onset of instability in a system combining uniform loss and local gain is a simple 
but new issue, for the comparison's sake it is relevant to briefiy consider it in a similar model, 
viz., the NLS equation: 

iut + ]^Duxx + \u\^u = -i'ju + (iFi - F2) S{x) ■ u, (15) 



where D is the spatial-dispersion coefficient, 7 and ri^2 having the same meaning as above. 
If the variable t in Eq. (|1^) is replaced by the propagation distance z, and x is reahzed 
as the transverse coordinate in a planar lossy waveguide, Eq. (|T5|) may be interpreted as 
describing spatial solitons in the above-mentioned case when the gain is applied along a 
narrow strip in the waveguide. In fact, the NLS model is a limit case of the BG system 
for small-amplitude solitons (see, e.g., Ref. |^); accordingly, Eq. (|T3|) is a small-amplitude 
limit of Eqs. (|^) and @. Nevertheless, it is pertinent to consider the NLS model parallel 
to the BG one, as it will help to understand the results for the BG system. 

A perturbation mode in the linearized equation (|l^) is sought for as [cf. Eqs. ) and 



u{x, t) = exp{—ixt — Kx) at x > 0, u{x, t) = A_ exp(— z^t + kx) at a; < 0, (16) 



with Re/t > 0. The substitution of Eq. (|T6D into the linearized version of Eq. (15) yields 



X 



(17) 



and the integration of Eq. (15) in an infinitesimal vicinity of x = leads to = A^, and 



Dk = r2 — iTi. 



(18) 



Note that the necessary condition Re k > and Eq. (|I^) show that, in fact, the perturbation 
mode ([TBD exists only in the case r2-D > 0. It is easy to understand the meaning of the 
latter condition: the inhomogeneity is attractive in this case, hence it can support the local 
mode. 

The substitution of Eq. (0) into Eq. (|17D yields a final result. 



x = - 



(TWTi 
2D 



+ ^7 



(19) 



As it follows from Eq. 

case. 



the instability-onset condition, Imx > 0, means, in the present 



Ti > (El 



(20) 



Thus, as well as in the case of the hot spot in BG, Eq. (^Of) demonstrates that the hot spot 
in the NLS model cannot give rise to the instability, unless it contains the imaginary part. 
Unlike the BG model, the additional condition r2-D > is also necessary for the instability. 
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C. An exact solution in the conservative model 



An exact solution for the pinned soliton is available for the conservative version of the full 
nonlinear model, with 7 = Fi = 0. In this case, it is easy to see that Eqs. (|^) and (|^) admit 
an invariant reduction, V{x) = —U*{x) (the asterisk stands for the complex conjugation), 
which leads to a single equation. 



i-^U + T2 cos6 + S{x) 
ax 



U + 3\U\^U-U* = 0. 



(21) 



As it follows from the integration of Eq. (^) around the point x = 0, the solution must 
satisfy a boundary condition 



U{x = +0) = U{x = -0) ■ exp(ir2) . 



(22) 



An exact soliton-like solution to Eq. (^l]) , supplemented by the condition (^21) , can be 
found, following the pattern of the exact solution for the ordinary gap soliton in the model 
with r2 = §: 



U{x) 



sin^ 



cosh 



[x + a sgn x) sin — ^9 



(23) 



where sgnx = ±1 for positive and negative x, and the real parameter a is determined by 
the relation 



tanh (a sin^) 



tan (12/2) 



tan(^/2) 

A corollary of the expressions (^Sf ) and (H), that will be used below, is 

,2 2 



|[/(x = o)r 



(cos To 



cos I 



(24) 



(25) 



From Eq. (|^) it follows that the solution exists not in the whole interval < ^ < tt, 
where the ordinary gap solitons are found, but in a region determined by the constraint that 
|tanh(asin6')| < 1, i.e., |tan(F2/2)| < tan (61/2), or 



iFol <9 <7i 



(26) 



(which implies that the solutions exist only if IF2I < vr). 

Although the exact solutions found above exist for either sign of F2, it is expected that 
only in the case F2 > they may be stable, as in this case the local inhomogeneity attracts 
the soliton (which is natural, as positive F2 corresponds to a local enhancement of the 
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refractive index, and a bright soliton is always attracted to an optically denser spot) [|T4 



In particular, in the case of small r2 the soliton may be regarded as a quasiparticle in an 
effective inhomogeneity-induced potential 

sin^ 6 . . 

= ^cosh (2esin^) +cos^ ' ^ ' 

where ^ is a displacement of the soliton's center from the point x = It is obvious that 

this potential indeed corresponds to the attraction and repulsion in the cases r2 > and 

T2 < 0. 

Note that the exact solution is a single-humped one, with a maximum at the point a; = 0, 
if r2 > 0; in the opposite case, the solution is a double-humped, with a local minimum at 
a; = and local maxima at x = ± |a|, as in this case Eq. (^) gives a < 0. In the limiting 
case 6 — \V2\ — > [see Eq. (pG])], Eq. (p5|) shows that \U{x = 0)| vanishes, i.e., the soliton 
pinned by the attractive inhomogeneity, with r2 > 0, reduces to zero, while the unstable 
two- humped state pinned by the repulsive inhomogeneity, with r2 < 0, goes over into a pair 
of two infinitely separated solitons with 6 = ~T2- 



D. The first-order approximation for the full model 



In the case 7 = Fi = 0, Eqs. (||) and (0) conserve the net energy. 



E 



\u{x)f + \v{x)\ 



dx . 



(28) 



(29) 



In the presence of the loss and gain, the exact evolution equation for the energy is 

^ = -27E + 2ri \\u{x)\^ + Hx)\^] u=o • 

at L J 

If the coefficients 7 and Fi are treated as small perturbations, the balance condition for the 
energy, dE/dt = 0, may select a particular solution, from the family of the exact solutions 
(p3|) of the conservative model, which remains, to the first approximation, a stationary 
pinned soliton in the full model. 

The balance condition following from Eq. (p9D demands 



7E = Fi \U{x = 0)f + \V{x 



(30) 



Substituting, in the first approximation, the unperturbed solution ( p3D and (|2^) into Eq. 
(0), and taking into regard the definition (pSf), the balance condition can be cast, after 



10 



some algebra, in a simple form: 



(31) 



cos r2 — COS 9 7 

As the pinned soliton may only be stable if r2 > 0, we consider this case. Note also that, 
according to Eq. (P^), we should constrain the consideration to the interval > T2, as 
otherwise the pinned solitons do not exist in the zero approximation (7 = Fi = 0). 

The pinned soliton selected by Eq. (|3ll) is expected to appear, with the increase of the 
stress parameter Ti/'j, as a result of some bifurcation. The inspection of Fig. 1, which 
displays {9 — F2) vs. F1/7, as found from Eq. (0), shows that the situation is qualitatively 
different in the cases F2 < 7r/2 and F2 > vr/2. 

In the former tangent (saddle-node) bifurcation [|1^ occurs at a minimum value 

(Fi/7)^j^ of the stress parameter at which Eq. ( pi]) has a physical solution for 9, and two 
solutions exist for F1/7 > (Xi/j)^^^. An additional analysis of Eq. ( |3lD demonstrates that, 
with the variation of F2, the value {'^i/'j)aim attains an absolute minimum, F1/7 = 1, at 
F2 = 7r/2. 

With the increase of F1/7, the lower solution branch that starts at the saddle-node bifur- 
cation point [see Fig. 1(a)] hits the limit point 6* = F2 [see Eq. (P^Dj, where it degenerates 
into the zero solution, according to Eq. (|25|). Equation (pT]) shows that this happens at 
the point F1/7 =1/ sinF2. On the other hand, it was shown above [see Eq. ([I^ ] that, 
precisely at the same point, the zero solution becomes unstable, in the limit of small Fi and 
7. According to the general stability-exchange principle |19|], the fact that the zero-solution 



branch gets unstable after its collision with another solution branch implies that the latter 
branch was already unstable. Therefore, we conclude that the branch which appears at the 
saddle-node bifurcation and ceases to exist hitting the zero solution, is an unstable saddle. 

The other (upper) branch generated by the saddle- node bifurcation [Fig. 1(a)] continues 
until it attains the maximum value 9 = 71 relevant to the physical solutions, which happens 
at 

£l=fEl) (32) 

This branch corresponds to the node-type solution which appears at the saddle-node bifur- 
cation point, therefore it has a chance to be stable. However, it may be unstable against 
perturbations that are not taken into regard by this elementary consideration (for instance, 
the possibility of a shift of the soliton's center from the point x = was not taken into 
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regard). The actual situation with the stabihty of pinned sohtons in the model including 
the loss and gain is rather complicated, see the next section. 

In the case r2 > 7r/2, the situation is different, as the saddle-node bifurcation is imag- 
inable in this case, occurring in the unphysical region 9 < T21 see Fig. 1(b). The only 
physical branch of the solutions appears at the point ri/7 =1/ sinr2, where it crosses the 
zero solution, lending it the instability, which the branch presumably had in the unphysical 
region. The stability-exchange principle which was already mentioned above suggests that, 
from the viewpoint of the present analysis, this branch becomes stable when it crosses into 
the physical region, 9 > T2- However, as well as the other branch considered above for the 
case r2 < vr/2, the present one may be subject to instabilities of other types. This branch 
ceases to be a physical one at the point (^). 

At the border between the two generic cases considered above, i.e., at r2 = vr/2, the 
saddle-node bifurcation occurs exactly at the point 9 = Tr/2, see Fig. 1(c). In this case, the 
destabilization of the zero solution happens at the same point. 

The situation is different in the case r2 = [see Fig. 1(d)], when the hot spot has no 
refractive- index-perturbation component, and Eq. ( pi]) takes the form 

2sin2 (9/2) ^ 'i ■ ^^^^ 

In this it was stressed above, the zero solution is never destabilized, in accordance 

with which the solution branches do not cross the axis 6* = in Fig. 1(d). The lower branch, 
which asymptotically approaches the ^ = axis, must be unstable (this is a generic feature 
in the case when the amplitude of the solution decreases with the increase of the stress 
parameter |]l9l), hence the upper branch may be stable within the framework of the present 
analysis. However, direct numerical simulations presented below demonstrate that, in the 
case r2 = 0, the pinned soliton is always unstable against the displacement of its center 
from the point x = 0. 

It may be relevant to compare these results with those that can be obtained for pinned 
solitons in the NLS model containing the loss and "hot spot", see Eq. (0). NLS solitons 
may only exist ii D > 0, therefore we adopt the normalization D = 1 in Eq. (plSj), and it 
makes sense to consider only the case when the inhomogeneity is attractive, i.e., r2 > 0, 
otherwise the pinned soliton has no chance to be stable. Note that, once we choose D > 
and r2 > 0, the zero solution may be unstable, according to the results presented above. 



12 



In the conservative limit, 7 = Fi = 0, the pinned NLS sohton is given by a commonly 
known solution, 

u = T^sech (rj {\x\ + a)) exp (^2^^^^) ' ^^^"^ 

"-^'■'(^^)' 

where rj is an intrinsic parameter of the solution family. Note that, as we assume r2 > 0, 
Eq. (pSD yields a > 0, hence the expression (0) has a single maximum at x = 0. 

If now 7 and r2 are introduced as small parameters, the energy-balance condition for this 
solution can be easily cast in the form rj = — r2 + 27/ri. In view of the relation (pSj), the 
actual solution exists in the region 77 > or, eventually, in the interval 

< Ti < 7/12 . (36) 

Comparing this result with the expression (^) that determines the instability threshold 
for the zero solution, and taking into regard that D = l,we conclude that the pinned soliton 
singled out by the balance condition disappears [crosses into the unphysical region, cf. Fig. 
1(a)], with the increase of the stress parameter, ri/7, exactly at the point where the zero 
solution loses its stability. According to the above-mentioned stability exchange principle, 
this implies that the zero solution inherits its instability from the soliton, hence the soliton 
solution, given by Eqs. (|3^) and (|35|) , is definitely unstable in all the region of its existence. 

This conclusion demonstrates that the above results for the pinned solitons in the BG 
model are nontrivial, as they give the pinned gap soliton a chance to be stable, which is not 
possible at all in the simpler NLS model. Actual stability of the pinned gap solitons will be 
studied below by means of direct simulations. 

III. NUMERICAL RESULTS FOR PINNED SOLITONS 
A. The approximation for the delta- function 

For the simulations, we have to adopt a numerical form of the 5-function in Eqs. (|l]),(0) 



and (^, (P). We use the same numerical scheme as in the recent work [Tj]. The scheme 
discretizes the coordinate x by 501 grid points Xj, j = —250, —1, 0, +1, ... + 250. As an 
approximation to the (5-function, the following function is defined on a set of 2A^ + 1 grid 
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points in the central part of the integration domain, located symmetrically around zero, 

J(.„_,„«,).^"= (37) 







elsewhere. 



The normalization factor A is defined so as to maintain the canonical normalization of the 
^-function, 6{x)dx = J2j ^ i^j) = 1; which yields 



A 



Ax cos I r^T : TT 

n=l 



2N +1 



(38) 



Ax being the spacing of the grid (in fact. Ax = 0.04). In most cases presented below, we 

(l + Vs) Ax] 



use = 2 [then Eq. (H) with Ax = 0.04 yields A 
makes the 5-function quite narrow indeed. 



7.726], which 



B. Stability of the pinned solitons in the conservative model 

Since exact stationary solutions to Eqs. (^ and (^ for the pinned soliton are available in 
the case 7 = Fi = 0, in the form of Eq. (p3|) supplemented by Eq. (p4|) , numerical test of their 
stability is straightforward. We simulated the stability by means of the split-step method 
applied to Eqs. (0) and @, employing the fast Fourier transform. The exact solution (^Sf) 
was taken as the initial configuration, and the corresponding value of 6 was varied. The 
values 6'in < IF2I, at which the exact solution does not exist [see Eq. (^6|)j were probed too. 
In this case, Eq. (|2^ ) yields an imaginary value of a, and the initial configuration was taken 



in the form of Eq. (23) with the imaginary a. Even though the latter configuration is not a 
stationary solution, it is still nonsingular and localized, so it can be used to launch the PDE 
simulation. 

As expected from what was mentioned above, in the case F2 < all the pinned states 
of the solitons are found to be unstable. Solitons are pushed away from the point x = 0, 
in accord with the expectation that the inhomogeneity is repulsive. It was also observed 
that, as IF2I increases, at F2 < —0.7 a small soliton is left behind around the point x = 
after the main pulse has separated from it; however, the residual soliton is also unstable and 
gradually decays into radiation. 

For positive F2, typical results regarding the stability of the pinned solitons are displayed 
in Fig. 2. A conclusion is that there is a single value ^stab ~ 7r/2 of the soliton parameter 
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6, such that if 6'in < 6 stab, the sohton decays into radiation, as is seen in Fig. 2(a). Sohtons 
with 6'in > ^^stab fclax into a stable one with 6 = Ogtah, see Fig. 2(b). Finally, a soliton with 
6'in = 6'stab directly gives rise to the stable soliton, see Fig. 2(c). 

The examples shown in Fig. 2 pertain to r2 = 0.4, and similar results were obtained for 
other values of Available computational power imposes a limitation on accuracy with 
which 6'stab can be identified. However, it was found that, for r2 = 0.1, the decrease of 
the soliton's amplitude, which is defined as \u{x = 0)|, is less than 1% after the evolution 
time T = 2007r, if 6^^ is taken from the interval (0.497r < < 0.527r), hence, in any case, 
6'stab (r2 = 0.1) belongs to the same interval. For a much larger value of the perturbation 
parameter, r2 = 1.1, the corresponding interval is 0.5l7r < < O.SSvr, hence 6'stab (r2 = 1-1) 
belongs to this region. Generally, 6'stab slightly increases with 

Figure 3 summarizes these results in the form of a plot in the (r2, 6'in) plane, which shows 
the regions where the initial soliton relaxes to the stable one or decays into radiation. In 
the region < \T2\, where the initial configurations are not true stationary solutions, this 
configuration decays into radiation immediately. 

These results, obtained for the conservative model with the local inhomogeneity of the 



refractive index, are very similar to those reported in Ref. ||T^ for the stability of the solitons 
pinned by an attractive inhomogeneity in the form of a local suppression of the Bragg grating. 
A noticeable common feature of the results is the existence of the single (up to the numerical 
accuracy available) value ^stab ~ 7r/2 of the parameter 6 which the established soliton may 
assume. In both conservative models (the ones considered here and in Ref. fl^), 6'in relaxes 
to 6'stab if 0in > 6'stab, and the soliton decays into radiation if 6'in < 6'stab, i-e., the soliton 
with 6 = Ostah may be called a semi-attractor. In fact, it strongly resembles semi-stable 
solitons, which are stable against small perturbations in the linear approximation, but may 
be unstable if terms quadratic in the perturbations are take into regard. Semi-stable solitons 



were recently studied in another context, as the so-called embedded solitons, see Ref. ||20 
and references therein. 

The fact that all the solitons with ^in > ^^stab, where ^stab is slightly larger than it/2, relax 
to the value 6 = Ostah, may be related to a known property of the ordinary solitons in the 
unperturbed BG model (7 = Fi = Fa = 0): they are unstable if ^ > 6^^^ ^ 1.011 ■ (7r/2) 
Thus, at least in the case when r2 is small, it is natural to expect that any pinned soliton 
with 9 > n/2 will relax, as a result of the instability, to a value close to 9^^\ What is less 
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obvious, is the decay of the sohtons with 6 < ^stab, and the fact that 6'stab so weakly depends 
on T2 (see Fig. 3). 



C. The pinned soliton in the lossy medium with the locahzed gain 

In direct simulations of the full model, which includes the loss and local gain, the exact 



solution (|23|) of the conservative version was again used as the starting point. The evolution 
of the solution was simulated at a fixed value of the loss parameter 7. The local gain Fi 
was varied in order to determine its value (s) at which the soliton settles down to a stable 
pinned soliton. 

Figure 4 shows the evolution of the soliton's amplitude, defined as \u{x = 0)|, vs. t, when 
the localized gain Fi is varied. The other parameters are fixed, so that 

T2 = 0.5, 7 = 0.0316, = O.Stt. (39) 

For a small value of Fi (Fi = 0.04208 in Fig. 4), which is insufficient to balance the loss, 
the soliton decays. For a slightly larger Fi = 0.04209, the soliton's amplitude grows, then it 
temporarily settles down (at the value 1.47 in Fig. 4), and, eventually, regular oscillations 
set in. A long simulation, up to t = 6007r (see Fig. 4) shows that the intrinsic vibrations of 
the soliton are completely stable. The waveforms |-u(x,t)| and |f(x, t)|, obtained at the end 
of the simulation for Fi = 0.04209, are shown in Figure 5(a). 

When the gain Fi is further increased, the initial growth of the soliton's amplitude is 
sharper; however, it is found that it again temporarily settles down to a nearly constant 
value close to the same level of 1.47 as above, which is followed by the onset of persistent 
oscillations. When Fi is still larger, the eventual oscillatory state becomes chaotic with large 
fluctuations. The corresponding waveforms of \u{x)\ and \v{x)\ at the end of the simulation 
{t = 300-7?) for Fi = 0.057 are shown in Fig. 5(b). It can be seen that conspicuous radiation 
tails are attached to the soliton, which implies a permanent energy leakage from it. This 
extra loss adds up to the direct dissipative loss, both being compensated by the localized 
gain. If Fi is too large, the radiation wave field outside the main pulse grows to such an 
extent that the resulting waveform cannot be regarded as a localized one. In fact, in this 
case it is observed that the main pulse separates from the point x = 0, drifts away, and dies 
down due to the loss. However, the strong localized gain generates a new "soliton" around 
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X = 0, which later drifts away again, this process repeating itself quasi-periodically. 

An important feature of these results is that a stable (even though it is vibrating) soliton 
is possible not at a single value of the gain, that exactly compensates the loss, but in a 
finite interval of values of the gain. The energy balance is maintained, in this case, through 
permanent emission of radiation by the soliton, which compensates the excessive gain. It is 
relevant to mention that a very similar mechanism, which gives rise to stable nonequilibrium 
solitons in an overpumped system of a different type (that, however, also originates in 
nonlinear optics - the so-called split-step model), was recently considered in detail in Ref. 



2^ . In that case too, the stability of the soliton is provided by the emission of radiation 



that balances the excess gain. 

A further insight in sustained intrinsic vibrations of the pinned soliton, and the transition 
from the regular oscillations to dynamical chaos, is provided by consideration of the spec- 
trum of the function \u{x = 0,)f:)|. In the estabhshed oscillatory regime, the spectrum was 
computed at several different values of Fi, while the other parameters were kept constant 



as per Eq. (|39D. Figure 6(a) shows the spectrum for Fi = 0.04209, which is the value barely 
enough to compensate the loss. It can be seen that the established oscillations are quasi- 
harmonic, with a single well-pronounced frequency 2.9 (in arbitrary units), and an additional 
tiny spectral component at the frequency ~ 2 (which is, apparently, incommensurate with 
the main one). 

Figure 6(b) shows the spectrum for Fi = 0.5633, which is similar to that in Fig. 6(a). The 
main frequency shifts down to a value about 2.8, with two other visible components found 
at the frequencies 0.8 and 1.4. Then, suddenly, at a slightly larger gain, Fi = 0.5634, many 
new conspicuous spectral components emerge, which is shown in Fig. 6(c), and corresponds 
to (apparently) chaotic intrinsic vibrations of the established soliton. The same behavior is 
observed at Fi = 0.5635. At Fi = 0.5636, the picture suffers another abrupt change [see Fig. 
6(d)]: the power spectrum reverts back to the simple three-frequency-component structure 
reminiscent of the situation at lower Fi, cf. Fig. 6(b). A transition from a chaotic behavior, 
(presumably) accounted for by a strange attractor, to a simple quasi-harmonic behavior is 
known in the theory of dynamical systems, where it is called a "boundary crisis" of the 
chaotic attractor [pl|]. 

The picture revealed by the simulations changes the third time at Fi = 0.5640, with 
reappearance of a many-component chaotic-like spectrum, similar to that in Fig. 6(c). The 
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chaotic behavior continues to higher values of Fi. Figure 6(e) shows the spectrum at Fi = 
0.05690, where its structure is not simply a mult i- component one, but continuous, which is 
characteristic for well-developed dynamical chaos. 

Another way to describe basic properties of the pinned solitons in the full model is to 
identify, for various initial values of ^in, a minimum value (Fi)^^;^^ of the gain which is 
necessary to overcome the loss. Figure 7 shows the evolution of the soliton's amplitudes as 
a function of time for 6'in = 0.2tt, O.Stt, and O.Qvr, the corresponding minimum values being 
(Fi)j^j^ = 0.0239, 0.0133, and 0.0131, while the other parameters are fixed, 7 = 0.01 and 
F2 = 0.5 [note that all these values of 6'in exceed F2, hence the corresponding exact solitons 
in the conservative model do exist, according to Eq. (|26|) ]. Thus, unlike the conservative 
model, in the full model, values of 6'in essentially smaller than tt/2 may give rise to a stable 
pinned soliton (with intrinsic vibrations). However, the smaller the difference 6'in — F2, the 



larger value of Fi is necessary, as, according to Eqs. (^) and (|25|), the rate at which the 
localized gain supplies energy to the soliton decreases ~ (^^in — F2) as 6'in — F2 0. On the 
other hand, analysis of the simulation results shows that the characteristics of the established 
soliton do not depend on the initial value of ^in which excited it, but solely on the values of 
7 and Fi 2, i.e., the established soliton is a genuine attractor. 

Then, effects caused by varying the loss parameter 7 were investigated. Because of 
the necessity to satisfy the energy balance condition, Fi needs to be changed to track 
the variation of 7. For each value of 7, we tried to find the minimum size of Fi that 
supports a stable soliton. Results of these numerical experiments, obtained for fixed 
F2 = 0.5 and ^in = 7r/2, and a set of values 7 = 0.000316, 0.001,0.00316, 0.01, 0.0316, and 
0.1, are displayed in Fig. 8. The respective minimum-gain values were found to be 
(Fi)^;^ = 0.000422, 0.00140, 0.00422, 0.0133, 0.04209, and 0.1327. It is interesting to note 
that, except for the second case, when the ratio {'^1)1^1^/'^ 1-40, in all the other ones the 
ratio takes values between 1.32 and 1.34. 

It is clearly seen from Fig. 8 that the amplitude of the established soliton monotonically 
increases with the growth of 7 (which is accompanied by the growth of the minimum gain 
necessary to support the soliton). It is also seen that it was never possible to produce a 
truly stationary soliton, but in some cases (7 = 0.000316, 0.00316, 0.01, and 0.1) it was 
possible to generate nearly stationary solitons with a small amplitude of intrinsic vibrations. 
In other cases (7 = 0.001 and 0.0316), varying Fi by steps as small as it was admitted by 
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the numerical scheme, it was not possible to adjust the gain so that to suppress the internal 
vibrations, i.e., the established soliton remained a breather, rather than anything close to a 
fixed-point state. 

In all the cases presented in Fig. 8, 7 and Fi are small enough to treat them as pertur- 
bations. Then, if the established soliton assumes a nearly stationary form, it is natural to 
expect that it must be close to the solution ( P5| ) found in the conservative model, with some 
value of 9, and this 9 must be related to 7 and Fi as per Eq. (|3l|) . 

It was checked that the quasi- stationary solitons, in those cases in when they were found, 
are indeed close to the wave form (^3]). The corresponding values of 9 were identified by 



means of the least-squared-error fit to the expression (23). Then, for thus found values of 9 
and given values of 7, the equilibrium values of the gain Fi were calculated as predicted by 
the analytical formula (|3T|) . Results of this are presented in Table 1. 

A noticeable fact obvious from Table 1 is that, in all the cases, the numerically found 
equilibrium value of the gain exceeds the analytically predicted one by 9 to 14 per cent. 
Because in all the C3iSGS, clS it was stressed above, the established solitons are not completely 
stationary, a natural conjecture is that the slightly vibrating soliton continuously emits 
energy at a low rate, and this extra energy loss makes it necessary to have the gain somewhat 
larger than that which compensates the direct dissipative loss as per Eq. (jHTj). 

As concerns the comparison of the full model with its conservative counterpart, we recall 
that, in the conservative model, the stable pinned soliton always assumes a single value of 
9 for given F2 (and this value very weakly depends on F2, always being slightly larger than 
7r/2, see Fig. 3). On the contrary to this, in the full model the quasi-stationary soliton may 
be stable in a range of the values of 9, as it is evident from Table 1. 

Finally, it has also been checked whether stable pinned solitons can be found when the 
"hot point" does not perturb the refractive index, i.e., F2 = 0. As a result, it has been 
concluded that any finite positive F2 (the smallest value tried was F2 = 0.01) may support 
a stable soliton in the pinned state, but if F2 = 0, the pulse set at x = finally drifts away, 
and then decays due to the loss. An explanation to this finding may be that all the solitons 
found in the model with loss and gain emit some radiation, see above, and asymmetric 
fiuctuations in the emission rate create a weak random force that drives the soliton away. 
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D. Soliton pinned between two repulsive inhomogeneities 

The soliton may be stably pinned not only by an attractive inhomogeneity, but also 
between two separated repulsive ones (in the present context, each one will represent a 
locally suppressed refractive index, corresponding to r2 < 0). The consideration of this 
configuration is interesting by itself, and it also may be used to design a soliton-based optical 
oscillator, in which the eigenfrequency is easily controlled by the choice of the separation 
between the two repulsive points. In particular, in the framework of the perturbation theory 
(for small |r2|), the soliton may be regarded as a quasiparticle in the external potential 

W^2(0 = -dj+w, + dj , (40) 

where the potential Wi{^) is given by Eq. (^) (with r2 < 0), and L is the separation 
between the two defects. 

We simulated the dynamics of this configuration in some detail, but only for the conser- 
vative case, 7 = Fi = 0. First of all, if L is smaller than the proper size of the soliton, it 
sees the pair of the inhomogeneities, in the first approximation, as a single repulsive cen- 
ter, hence stable bound states are not possible. Within the framework of the perturbation 
theory, Eqs. ( ^OD and ( pTj) make it possible to predict a critical value, (A^)^^., at which 
a stable equilibrium appears for the first time at x = 0. The corresponding expression is 
cumbersome, but it easy to verify that (A,^)^^, monotonically decreases, with 9 varying from 
to 7r/2, from (A^cr = oo to the minimum value (A^^ = In (v^ + l) ^ 0.88. 

Direct simulations at finite F2 demonstrate [see an example in Fig. 9] that, in the case of 
relatively small L, when the pinned state of the soliton is unstable, the dynamical evolution 
does not trivially reduce to pushing the soliton aside; instead, a generic outcome is splitting 
of the soliton in two, which is accompanied by a spontaneous symmetry breaking (in some 
cases, for instance if 9\a = O.Tvr, the other parameters being the same as in Fig. 9, the 
instability develops so slowly that it was not possible to identify the outcome). 

With the increase of L, stabilization of the soliton trapped between the repulsive inhomo- 
geneities becomes possible. The trapped states seem most stable around the value L = 3.84, 
see an example in Fig. 10. In this case, systematic simulations reveal a feature which 
strongly resembles the one reported above for the single attractive inhomogeneity in the 
conservative model: an established trapped state is stable for a single (up to the accuracy 
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of numerical simulations) value of 9, which is very close to 7r/2; if > 7r/2, the soliton 
sheds off some radiation and eventually relaxes to the said single value of 6 (see Fig. 10), 
while if ^^in < 7r/2, the soliton gradually decays into radiation. Thus, the single- valuedness 
of the stable soliton in the conservative model appears to be a generic property. For still 
larger values of L, the pinned state is less robust; in particular, a soliton with 9 > n /2 
may split, instead of relaxing to the stable one with 6 ^ tt/2. 

IV. COLLISION OF A MOVING SOLITON WITH THE LOCALIZED GAIN 

Once the existence of stable pinned soliton has been established, the next natural step is 
to consider a possibility of capturing a free moving soliton by the "hot spot" . To this end, 
the soliton was first generated far from the spot by means of the Newton-Raphson method, 
as a stationary solution in the reference frame moving at some velocity c; a range of the 
velocities < c < 0.7 was thus investigated. Then, the collision was considered, running 
direct simulations of Eqs. (0) and (|^). 

First, the collision experiment was performed in the conservative model, with Fi = 7 = 0. 
The parameter plane (c, F2) was explored with c taking values 0.1, 0.2, . . . , 0.7, and F2 taking 
values 0.1, 0.2, . . . , 0.9, while was kept constant at 0.77r. 

If the inhomogeneity is weak, the moving soliton passes through it, see an example in 
Fig. 11(a). When the inhomogeneity strength is larger, r2 0.5, a part of the soliton still 
passes through it, but another part of the soliton's energy is captured by the local defect to 
form a pinned soliton, an example of which is shown in Fig. 11(b). Some radiation bouncing 
in the backward direction can also be observed when r2 is large, or when the incident soliton 
is fast. Naturally, more energy is trapped by the defect if r2 is larger [Fig. 11(c)], and less 
energy is trapped if the soliton is faster. However, the value r2 ~ 0.5, at which the trapping 
begins, only weakly depends on the soliton's velocity c. Figure 12 summarizes these results, 
showing a border in the (c, r2) plane between the passage and partial-capture regions. 

Next, we consider the collisions in the full model, with 7 = 0.01 and Fi = 0.015. Results 
reported in the previous section show that a stable pinned soliton exists at these values of 
the loss and gain (the collisions were simulated only for small values of 7, as otherwise the 
soliton will be strongly attenuated still before the collision). The initial value of the soliton's 
parameter was again = O.Tn. 
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The results obtained for the full model are not very different from those for the conser- 
vative one. When the gain r2 is small, the soliton passes through, and if r2 is larger, a part 
of the energy is trapped to form a pinned soliton. A difference from the conservative model 
is that the value of r2 at which the inhomogeneity starts to capture a part of the soliton's 
energy in the conservative model is approximately independent of its velocity: r2 ~ 0.5 if 
c > 0.1, while in the lossy model, this value of r2 increases with c, as is seen in Fig. 12. 

Another representative set of numerical data can be displayed for a fixed value of the 
soliton's velocity, c = 0.1, while the parameter of the incident soliton takes values 
O.lvr, 0.27r, . . . , O.Qvr, and V2 = 0.1, 0.2, . . . , 0.9. These results are presented here only 
for the conservative model. 

Simulations show that, if the moving soliton is heavy (large 6[^^) or the inhomogeneity 
is weak, the soliton passes it. Heavier solitons can pass through a stronger defect. If the 
inhomogeneity is strong (r2 is large), the soliton gets trapped, which is always accompanied 
by emission of radiation in both the forward and backward directions, and the radiation 
can further self-trap into secondary solitons. At small r2, little energy is scattered away in 
either direction. If r2 is larger, more energy is scattered forward, and when r2 is still larger 
(r2 ~ 0.9), more energy is scattered backwards, cf. Fig. 11(c). Figure 13 summarizes the 
results obtained for the interaction of the moving soliton and the localized attractive defect 
in the conservative model for the fixed velocity, c = 0.1 

It is natural to compare the results obtained for the conservative model with those re- 



ported in Ref. |1J] for the interaction of the moving gap soliton with an attractive inho- 
mogeneity in the form of a local suppression of the Bragg grating. In that case, when the 
soliton was heavy (large 6), the interaction effectively reverted from attraction to repulsion, 
so that the incident soliton could bounce back. In the present model, this unusual behavior 
has never been observed. 



V. CONCLUSION 

In this work, we have introduced a model of a lossy nonlinear fiber grating with a "hot 
spot" combining the localized gain and attractive inhomogeneity of the refractive index; the 
spot can be created by means of doping a short segment of the fiber. In the absence of the 
loss and gain, a family of exact solutions for pinned solitons was found. In the full model 
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including loss and gain, the instability threshold for the zero solution was found; it was 
concluded that the instability is not possible without the presence of the imaginary part of 
the local gain, i.e., a localized perturbation of the refractive index. Further, for small values 
of the loss and gain, it was predicted what soliton is selected by the energy-balance condition. 
Parallel to this, it was shown that, in the simpler model based on the NLS equation, the 
pinned soliton can never be stable in the presence of the loss and local gain. 

In direct simulations, we have found that a single pinned soliton is stable in the conser- 
vative fiber-grating model. It is a semi-attractor: solitons with a larger energy relax to it via 
emission of radiation, while the ones with smaller energy completely decay into radiation. 
The same conclusion is obtained for solitons trapped between two repulsive inhomogeneities. 
In the full model with the loss and gain, all the stable pinned pulses demonstrate persis- 
tent internal vibrations and emission of radiation. Sometimes, they are almost stationary 
solitons, and in these cases the prediction based on the energy balance underestimates the 
necessary gain by 9% to 14%, which is explained by the extra radiation loss. If the loss and 
gain increase, the intrinsic vibrations become chaotic. 

Collisions of free moving solitons with the "hot spot" were simulated too. The passage 
and capture regimes were identified for the solitons in the conservative and full models; the 
capture is only partial, which actually implies splitting of the soliton. If was also found 
that, if a large part of the soliton's energy is radiated away, it may self-trap into secondary 
solitons. 
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Tables 



7 


9 


("^1 )num 


(-'^l)anal 


(^l)num (^l)anal 

(ri)„n„i 


0.000316 


O.Stt 


0.000422 


0.000386 


0.0944 


0.00316 


0.5957r 


0.0042 


0.00369 


0.1373 


0.01 


0.6087r 


0.01333 


0.01165 


0.1442 


0.1 


0.8267r 


0.1327 


0.121 


0.0967 



Table 1. Values of the loss parameter 7 at which quasi-stationary stable pinned soli- 
tons were found by the adjustment of the gain Fi (see Fig. 8), while the refractive- index 



perturbation is fixed, F2 = 0.5. The corresponding values of the gain, (Fi 



'num' 



are also 



included, together with the values of the soliton parameter 6 which provide for the best fit 
of the quasi- stationary solitons to the analytical solution (|^). The values (Fi)^^^^^^^ are those 



predicted, for given 7 and 6, by the energy-balance equation (0), which does not take the 
radiation loss into account. 
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Figure captions 

Fig. 1. Analytically predicted solution branches for the pinned soliton in the case of 
weak loss and gain. Shown is — r2 vs. the stress parameter ri/7. (a) An example of the 
case r2 < 7r/2 is displayed with r2 = vr/4; (b) an example of the case r2 > 7r/2 is displayed 
with r2 = 3'7r/4; (c) r2 = vr/2; (d) r2 = 0. In the last case, nontrivial solutions appear 
at the point 9 = 0.74427r, ri/7 = 1.3801, and at large values of ri/7 the continuous curve 
asymptotically approaches the horizontal axis. In all the panels, the dashed lines show a 
formal continuation of the solutions in the unphysical regions, 6 < and 9 > tt. In the 
panels (a), (b), and (c), the trivial solution, 6' — r2 = 0, is shown by the solid line where it 
is stable; in the case corresponding to the panel (d), all the axis = corresponds to the 
stable trivial solution. (Note that all quantities plotted are dimensionless.) 

Fig. 2. Typical results illustrating the stability and instability of pinned solitons in the 
case 7 = Fi = 0, and F2 = 0.4. Each panel shows the evolution of \u{x, t)\, starting with the 
exact pinned-soliton configuration (the evolution of |f(x, t)| is quite similar.) (a) = OAir 
is smaller than 6'stab: the soliton decays into radiation, (b) Oi^^ = O.Svr > ^^stab^ the initial 
soliton transforms itself into a stable one, shedding off excess energy in the form of radiation, 
(c) din = 0-57T ^ ^stab- Direct appearance of the stable soliton. 

Fig. 3. A summary of results obtained for the stability of pinned solitons, plotted in 
the (F2,6'in) plane, in the conservative model with 7 = Fi = 0. In the upper region, where 
> 6'stab, initial solitons evolve into the stable one, shedding off extra energy. In the lower 
region, where 9 < 6'stab, solitons completely decay into radiation. Beneath the lower solid 
line, which borders the region where IF2I < 6 < n, see Eq. (|26|) , stationary solutions for the 
pinned solitons do not exist. Accordingly, an initial pulse taken as a formal "soliton", with 
an imaginary root of Eq. (^) substituted for a (see the text), is immediatley destroyed. 

Fig. 4. Evolution of the amplitude of the pinned soliton in the full model with loss 
and gain, in the case with F2 = 0.5, 7 = 0.0316 and 6-,^ = O.Svr. If Fi = 0.04208, the 
gain is insufficient to balance the loss, and the sohton decays. When Fi = 0.04209, the 
soliton initially grows, and its intrinsic vibrations develop. When Fi takes a slightly larger 
value, 0.04215, the initial growth of the amplitude is steeper, which again results in the 
establishment of regular intrinsic vibrations (in this case, the oscillations are very similar to 
those supported by Fi = 0.04209). When Fi is essentially larger, for instance, Fi = 0.057, 
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interal vibrations of the pinned soliton become chaotic (the latter case shown by the dotted 
hne) . 

Fig. 5. The profiles of \u{x,t)\ (sohd lines) and \v{x,t)\ (dashed lines) at the end of the 
simulation (narrow peaks placed at a; = mark the "hot spot"). The values of 7, and 6[,^ 
are the same as in Fig. 4. (a) Fi = 0.04209 is barely enough to support the soliton against 
the loss. In this case, the soliton emits radiation at a low rate, (b) Fi = 0.057. The soliton 
emits radiation at a high rate. 

Fig. 6. The frequency spectrum of the time-dependent amplitude \u{x — 0,t)\ of the 
pinned soliton at different values of the local gain Fi, after persistent vibrations set it: (a) 
Fi = 0.04209, (b) Fi = 0.05633, (c) Fi = 0.05634, (d) Fi = 0.5636, and (e) Fi = 0.05690. 
In all the cases, 7 = 0.0316, F2 = 0.5, and = O.bn. 

Fig. 7. The soliton's amplitude \u{x = 0,t)\ vs. t for three different values of the initial 
amplitude ^ in. In each case, the value of the gain Fi is chosen as the minimum one which 
can support the establishment of a soliton. Other parameters are fixed: F2 = 0.5, and 
7 = 0.01. 

Fig. 8. The amplitude of the soliton, \u{x = 0,t)\, vs. t for different values of the loss 
parameter 7. Each time, the value of the gain Fi is chosen as the smallest one which leads 
to the establishment of the soliton. Other parameters are fixed: F2 = 0.5, and — O.Btt. 

Fig. 9. The interaction of a soliton with a pair of repulsive points (F2 = —0.7), with 
a relatively small separation between them, L — 1.84, in the conservative model, (a) The 
initial configuration, with — 7r/2. (b) The result of the interaction: sphtting of the sohton 
into two pulses, which is accompanied by a spontaneous symmetry breaking. The solid lines 
show and the dashed lines show \v\. Note that, in the initial configuration, \u\ and \v\ 
are indiscernible. 

Fig. 10. Stable soliton captured between the repulsive points (F2 = —0.5) with the 
separation L — 3.84 between them, in the conservative model. Shown is the evolution of the 
field \u{x, t)\. The initial value of the soliton parameter is = 0.77r. 

Fig. 11. CoUision of a moving soliton, with fixed values 9 = O.Tn and c = 0.4, and the 
inhomogeneity in the conservative model (the inhomogencity is shown by a narrow peak 
which, for an unessential reason, is shifted from the point x = 0). The lower and upper 
panels show, respectively, the evolution of the field |K(a;, t)| in terms of the contour plots, 
and the waveforms \u{x)\ and \v{x)\ (solid and dashed lines) at the end of the simulation 
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(note that the u- and I'-components are asymmetric in the moving sohtons). (a) The sohton 
passes through a weak defect with r2 = 0.2. (b) A stronger defect, with r2 = 0.5, captures 
a part of the energy of the passing sohton, to form a small-amphtude pinned one. Another 
small part of the energy bounces back in the form of radiation, (c) If the defect is still 
stronger, r2 = 0.9, the shares of the trapped and bounced energy are larger. 

Fig. 12. Borders in the parametric plane (c, between regions in which the moving 
soliton with fixed 9 — O.Ttt passes the defect or gets partially trapped. The solid line is the 
border in the conservative model, with 7 = Fi = 0. The dashed hue is the border in the full 
model with 7 = 0.01 and Fi = 0.015. 

Fig. 13. Regions in the parametric plane {6, F2) of the conservative model in which the 
moving soliton with fixed c = 0.1 passes the defect or gets partially trapped. 
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